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Abstract 

The superconducting properties of a recently proposed phenomenological model for a weakly 
doped antiferromagnet are analyzed, taking into account fluctuations of the phase of the order 
parameter. In this model, we assume that the doped charge carriers can't move out of the anti- 
ferromagnetic sublattice they were introduced. This case corresponds to the free carrier spectra 
with the maximum at k = (±7r/2, ±7t/2), as it was observed in ARPES experiments in some of 
the cuprates in the insulating state The doping dependence of the superconducting gap and 
the temperature-carrier density phase diagram of the model are studied in the case of the d x 2_ y 2 
pairing symmetry and different values of the effective coupling. A possible relevance of the results 
to the experiments on high-temperature superconductors is discussed. 
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I. INTRODUCTION 



Theoretical description of high-temperature superconductors (HTSCs) remains one of 
the most important unsolved problems of modern condensed matter theory. However, de- 
spite the lack of complete understanding of this phenomena, some facts about the cuprate 
materials are almost generally accepted by scientific community. First of all, it is believed 
that superconductivity in this materials is mainly generated in CuC>2 layers, since even sam- 
ples with the layer spacing of ~ 10A demonstrate very high critical temperatures. Second, 
these materials in undoped regime are antiferromagnetic insulators, which transform into 
superconductors when the doping is larger than some critical value. Experiments indicate 
the presence of strong antiferromagnetic correlations in the superconducting regime 

as. 

Third, the symmetry of the superconducting order parameter is of a c/-wave type in many 
of the compounds 0[ |^ . It is believed by many researches that these facts contain sufficient 
information to understand the phenomena of HTSC. It is almost generally accepted that 
this phenomena can be described in the framework of a simple two-dimensional model of 
strongly correlated electrons (for a recent over-review, see jf]]). The antiferromagnetic cor- 
relations produce a BCS-type spin-wave pairing between electrons. One of the most often 
studied model of HTSCs is the t — J- model |3, S] (see, for example, 0, 10, 11, 12. llH ll^|). 
It was shown that the solution of this model within different approximations demonstrates 
c?-wave superconductivity in the underdoped regime. This and many other models, however, 
are rather oversimplified, and they do not take into account some properties of the cuprates. 
In particular, it is known from ARPES experiments that the free carrier spectra of some 
cuprates like S^CuC^C^ in the insulating phase has its maxima not at the Brillouin zone 
edge at momenta k = (±7r, ±7t), but at points k = (±7r/2, ±vt/2) [lj]. 

This fact suggests that the free carrier spectra must be e(k) ~ —4t nn cosk x cosk y — 
cos 2k x + cos2ky), i.e. the doped carriers move only within one sublattice where they 
were born. This situation is similar to the case of collinear antiferromagnetic dielectrics with 
the Neel ground state, where electronic and spin excitations can't move out of one sublattice 
to another. It is important that the excitations must live on non-equivalent sublattices, i.e. 
on the copper sublattices, which form an antiferro magnet. It is often assumed for several 
reasons (see below), that doped holes in cuprates occupy sites in the oxygen subsystem of 
the Cu02 planes However, it is known that the n.n.n. hopping parameter for this 
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sublattice is much smaller than the n.n. hopping parameter, which excludes the possibility 
to explain the spectra experiments [l] in this case. One can try to overcome this difficulty 
by introducing the Zhang- Rice singlet, which consists of a doped oxygen hole coupled to one 
of the closest copper ions. The hole chooses the ion in such a way that the energy of the 
singlet is minimal. Such states can move only within their magnetic sublattices. However, 
these states are unstable due to several reasons. First of all, the "up" and "down" spins of 
the copper ions enter in the linear combination of the Zhang-Rice singlet in an equivalent 
way, so a hole with the opposite spin projection will try to form the same singlet with 
the neighboring localized ion spin from another sublattice. Also, the oxygen ions occupy 
states which are symmetrical with respect to the neighbor copper ions. Therefore, the total 
exchange field which acts on a hole on the oxygen site is compensated, and the hole becomes 
frustrated with respect to the choice of the axis of the spin quantization. The proper axis of 
quantization can be established, but in this case the hole hybridized with the copper states 
will move in the Cu02 plane not feeling the magnetic ordering. This also doesn't allow one 
to explain the spectra experiments. 

These difficulties can be avoided by assuming that the doped holes sit on the copper sites 
[r|. In fact, in the undoped state the Cu0 2 planes consist of ions Cu 2+ and 2 ~. The 
doping leads to the valence change of some ions, which results in appearance of movable 
charges in the system. The assumption that the holes occupy the copper ions corresponds 
to the transitions Cu 2+ — > Cu 3+ , which are different from usually considered transitions 
2 ~ — > O in the case of the holes which occupy the oxygen sites. The reason for the last 
assumption follows from the fact that it is believed by many researches that the configuration 
Cu 3+ 2 ~ has higher energy comparing to Cu 2+ 0~. However, the configuration Cu 3+ 2 ~ can 
be preferable if one takes into account the Coulomb attraction Vc inside the configurations. 
In fact, since 

Y Cu»+o»- ^ 3V g« 2+ o- ~ 6Vg u+ °", and every copper ion is surrounded by four 
oxygen ions, at the same time as every oxygen is surrounded by only two copper ions, one can 
find that the corresponding Coulomb attraction energy difference for two configurations is 
of order 20Vq u+o . This difference in the Coulomb interaction energies can provide stability 
of the Cu 3+ 2_ configuration comparing to other configurations, i.e. a possibility for the 
doped holes to move within the copper sublattices. This picture corresponds qualitatively 
to the t — J-model case. However, it was usually assumed that the free carrier dispersion 
relation in this model corresponds to the nearest neighbor hopping, which does not allow 
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one to obtain the experimental spectra with the minimum at k = (±7r/2, ±7r/2). 

A simple phenomenological Hamiltonian, which corresponds to charge carriers that move 
within fixed sublattices, was proposed in order to describe superconductivity in cuprates at 
low carrier densities in |l6j. The effective attraction between the doped electrons on different 
sublattices was assumed to be equal to the antiferromagnetic coupling J between nearest 
site spins. In fact, the carrier doping in the antiferromagnet leads to a breakdown of the 
antiferromagnetic coupling J between nearest site spins. This results in an increasing of the 
energy of the system. This energy increasing is minimal when two empty sites are the near- 
est neighbors, because the minimal number of the antiferromagnetic spin-spin (exchange) 
interaction bonds is broken in this case. This phenomenological attraction was introduced 
for the first time by Trugman in Ref.Q]. As it was shown in Ref. [3], the order parameter, 
which corresponds to the pairs produced by such an interaction, has a <i-wave symmetry. 

Below, the doping dependence of the c?-wave superconducting gap in the model at different 
values of coupling and T = is studied. In addition, we analyze the temperature-doping 
phase diagram of the model. It is shown that the width of the pseudogap temperature 
region between the mean-field critical temperature T C MF and the critical temperature of 
the Berezinskii-Kosterlitz-Thouless (BKT) transition T c (= Tbkt) strongly depends on the 
carrier concentration. We make a qualitative comparison of the results with the experimental 
results on some HTSCs. 



II. MODEL 



The Hamiltonian of non-interacting doped d-hole carriers in HTSCs can be approximated 

by 

H d = £ d Yl Yl d ^Jna n - \ Y Y W ff nkm>dL n rf m«x m , (1) 



2 

n cr n n,m (T n ,&-n 



where 4u„(^nff„) is the creation (annihilation) operator of the electron on site n with spin 
cr n , Ed is the electron on-site energy, t nm is the hopping parameter and (a n \a m ) is spin-spin 



correlation function calculated in the ion system of coordinates (see |l6l|). This Hamiltonian 
can be transformed to the following form in terms of the Hubbard operators in the laboratory 
system of spin coordinates: 

H d = H^ + Hj»+H®, (2) 



where 



H coh = (e d - /,) E cos Q " FM( 2 n " m) ^ 1/2 ^ /2 ' 2 (3) 

n n,m 

is a part of the Hamiltonian which describes motion of free holes in an antiferromagnetically 
ordered medium, Qafm & = (i 71 ") i 71 "); a is the lattice constant of a square lattice, 

= 4 £ sin ^l 11 "^ W 1/8 ^ Aa 5- - (4) 
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n,m 



and 

fl2 = ~ E «- » s Q " FA ' ( 2 n " m) ^' 1/2 ^ /2 ' 2 g^ m (5) 

n,m 

describe non-coherent inter-ion hole transitions with one and two spin excitations, corre- 
spondingly. S£ and S m are spin creation and annihilation operators. It is important that 
expressions ©-© are written in terms of the Hubbard operators, which directly take into 
account the antiferromagnetic ordering in the system. The ion spin projections in both mag- 
netic sublattices are equal to 1/2 in the ground state of the crystal (we use local systems 
of coordinates for each sublattice). The spin conservation is also taken into account, which 
results in the fact that the electrons can move only on the magnetic sublattice on which 
they were born (for details, see 

As it was mentioned in the Introduction, a simple effective attraction between the doped 

n 

electrons on different sublattices can be introduced [18]. The carrier doping in the anti- 
ferromagnet leads to an increasing of the energy of the system, since it breaks the antifer- 
romagnetic coupling J between nearest site spins. The energy increasing is minimal when 
two doped particles occupy the nearest neighbor sites, since the minimal number of the 
antiferromagnetic spin-spin bonds is broken in this case. Therefore, the doping leads to an 
effective attraction between carriers on different sublattices: 

Hattr = -J 22 Xn' 2 Xn+p- (6) 
n,p=a,b 

The total Hamiltonian of the system is H = + H attr . For simplicity, we neglect terms 
(J1J and (0), which do not contribute to a significant increasing of the superconducting 
critical temperature. The first term corresponds to a BCS-like interaction, which is small 
at low carrier densities, or Fermi momenta of doped holes kp. It is proportional to the 



energy of spin waves ^afm^f), which is much smaller than the exchange energy J. The 
term Eq. (jSJ), which corresponds to two-magnon attraction, also doesn't contribute to a 
significant increasing of the critical temperature in the <i-wave pairing channel 17 1. 

Therefore, a simplified version of the Hamiltonian for the doped the antiferromagnet can 
be written as: 

H = (e d - £ X? - \ £ U cos QAFA/( 2 n " m) X^/^ - J £ X^ p . (7) 

n n,m n,p=a,b 

In the case of an antiferromagnet on a square lattice with two sublattices, the free particle 
energy spectrum, which corresponds to the first two terms in Eq. (JJJ), is: 

e(k) = Ed — 4£ 2 cos k x cos k y — 2t 3 (cos 2k x + cos 2k y ) — /x, (8) 

where /x is the chemical potential, and i 2 and £3 are the next nearest and next next nearest 
neighbor hopping parameters, correspondingly. We use the units where the lattice constant 
is equal to one: a = 1. In order to have e(k) = at k = one can choose Ed = 4t 2 + 4£ 3 . 
The chemical potential is connected with the free (or doped) particle number in the system 
by the following natural relation: 

* = J2(X 2 n \ (9) 

n 

where the sum goes over two sublattices. The Hamiltonian Eq. (JJJ) has a simpler structure, 
comparing to the t-J-model Hamiltonian, yet it can describe some of the main physical 
properties of underdoped cuprates. 

III. ZERO-TEMPERATURE PROPERTIES 

To study superconducting properties of the system described by the Hamiltonian Eq. (JJ|), 
it is convenient to introduce generalized Nambu-Hubbard hole operators: 

*-(*) = 1^2^ ] ' ^ = W^)'^' 1 ^)) , (10) 

where n are lattice sites and t is time. In this case the time ordered Green function 
G nm (t,t') = -i(T(^ n (t)^ m (t')) is 

G {f t>) = _ % . (TX 2 n 1/2 (t)X^ 2 (t')} (TX 2 n 1/2 (t)X 2 m 1/2 (t>)} 

' 1 (TX^^X^it')) (TXi /2 ' 2 (t)X 2 m 1/2 (t')) 1 ' 



(13) 



The Green function Eq. flllj) satisfies the following equation of motion: 

*^G nm (t, = 6{t - t')S a J + (T[* n (t), fl]*L(0> 5 (12) 
where, as it was mentioned above, if is defined by (J7J), and 

(X^ /2 ' 1/2 (t)+X^( t )) 

<x n /2 - 1/2 + Xf) 

This equality can be derived by using the commutation relations for the Hubbard operators. 

In the generalized mean-field theory approximation, the last term in Eq. (jl2[) can be 
linearized in the following way (see, for example Ref. [20]): 

(T[^ n ,H]¥j(uj) -J2^yG lm (u;), (14) 

l 

where 

JSU = <{[*», ^.*!a}> (15) 
is the energy matrix. The nonlinear (dynamical) corrections to the self-energy in Eq. (|14|) 
can be taken into account [12]. We assume that the generalized mean-field approximation 
Eq. ((Hj) is good enough in the case of low carrier densities, when the free quasi-particle 
excitations correspond to the field represented by the Hubbard operators Eq. ijlH]) . 

The expression for the energy matrix (|15|) can be found by solving the Heisenberg equa- 
tions of motion for the X-operators. In terms of the energy matrix Eq. (|15|). the Green 
function can be written as: 

G nm (u) = *_ . (16) 

To find the explicit expression for the Green function Eq. ([16)1 . let us write down the equa- 
tions of motion for the Hubbard operators: 

d_ 

at' 



-\ £ *m cos ^"^ (ir + X?)X}'V - 2 J £ X^X^ 2 , 
i p 

(17) 



i»^ 1/2 (*)= - (^-^)^' 1/2 
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dt 

+ \ £ *in cos ^(l-n) Xi2 , 1/2( ^ Al/2 + x , 2) + 2 7 ^ ^, 1/2 ^ 

1 p 

(18) 



Substitution of the expressions Eqs. (JTTj) and (|T%|) instead of the anti-commutators 
(idX/dt = [X, H]) into Eq. (fTB^) and evaluation of the anti-commutators give the follow- 
ing expression for the energy matrix: 

^nm E nm T z -\- A nm T.j., (19) 

where f x and f z are the Pauli matrices and 

P X ( r- ,A/vl/2,l/2 i Y-2,2\ | ^ / r,^ Q-A-FmO- — n ) / v-2,1/2 y l/2,2\ 
Km= - nm (£: d - fJ,)[X n ' '< +A n ' ) + -d nm ^ln COS (Aj A n ' ) 

1 

+ 25 nm £ J nl (Ay 2 'V2 x 2,2 ) + l tmn cos Q^M(m-n) ( ^ Al/2(x , Al/2) + ^ 

(20) 

A ^ X ™o QaFmO — n) 2,1/2 v 2,l/2\ 7 / v2,l/2y2,l/2\ / 1 \ 

A nm = --d nm 2_^iinCos (Aj ' X n 1 ) - 2J nm (A n - 7 7 ). (21) 

l 

are the renormalized energy and the superconducting gap matrices. Despite their compli- 
cated formal structure, it is possible to show that the terms in Eq. (|2Tlj) and (|21|) have a very 
simple physical interpretation. In particular, the second and the third terms in Eq. (|2(J|) . 
proportional to 5, lead to a renormalization of the chemical potential fi — > y! = fi + 5y. The 
average (X]i 2 ' l ^ 2 {Xn 2,1 ^ 2 ) + X 2,2 )) multiplied by the hopping operator t mn in the last term 
of Eq. (jlZUj) leads to renormalization of the quasiparticle bandwidth in the limit of low dop- 
ing;. However, one can put (X^ 2 ' 1/2 (Xi /2 ' 1/2 ) + Xl' 2 )) ~ 1 in this limit, since in the limit of 
low carrier concentrations xl/ 2,1 ^ 2 ) + X^ 2 ~ 1, and the renormalization of the quasiparticle 
band is not strong. Therefore, the energy function Eq. (j2U|) can be approximated by the 
free energy spectra expression Eq. (jHj) multiplied by —1 in the momentum space. It is also 
assumed that the chemical potential in Eq. (|SJ) is renormalized. 

The expression for the gap function Eq. ()21|) can be also simplified. In fact, as it was 
shown in Ref. the attraction (jUJ) favours a superconducting pairing with the <i-wave 
symmetry of the order parameter. Therefore, we assume that the strongest pairings in the 
system takes place in the <i-wave pairing case, and neglect the first term in Eq. (j'21|) . which 
doesn't contribute to the pairing in the d-wave channel. Therefore, the gap function can be 
approximated in the following way in the momentum representation: 

A(k) = - 2 £ J(k-q)(A 2 f Xji /2 ); (22) 
q 



S 



where we introduced the nearest neighbor attraction kernel J(k) = 
2 J [cos(k x — q x ) + cos(/c y — q y )]. Thus, Green function (fTBj) has the following form in 
the momentum space: 

GV> k ) = 7TT^ — 777^ ( 23 ) 

where we assumed that in the limit of low carrier densities the normalization matrix is 
approximately equal to the unit matrix J ~ 1. 

To find the unknown gap function A(k) and the renormalized chemical potential //, 
one can write down and solve the system of equations for these functions by using the 
fluctuation-dissipation theorem: 

. 1 f°° , duj%G AB (Lo) , lS 



The zero-temperature equations which connect A(k) and \l with the parameters t?,, £3, J 
and 5 follow from Eqs. ((22) and ©: 

A(k) = -2 £ J(k - q) =4==7 r (25) 

^ = v[l+ - £(k) 1. (26) 
^ 2 ( k ) + A2 ( k )J 

Since we consider the case, when the pairing in the system takes place in the ci-wave 
channel, we put A(k) = A^7^(k), where 7<z(k) = (cosfc^ — cosk y ), in Eqs. and In 
order to these equations in this case, one must extract the d-wave piece form the interaction 
kernel and approximate it by this function: J(k — q) — > 2 J7^(k)7rf(q). In this case, Eq. ()25|) 
for the superconducting order parameter acquires a rather simple form 

l=4J^ 7 , 2 (q) V 2 , .]. 22( ■ (27) 

q V£ 2 (q) + A^(q) 

The solution of the set of Eqs. (j2T?|) and (|27j) at different values of interaction and 
hopping t 3 is presented in Fig.l. As it follows from this Figure, the gap is not very sensitive 
to the values of the next-nearest neighbor hopping parameter, but it strongly depends on the 
interaction potential. Superconductivity is suppressed when the carrier density is smaller 
than some critical value. This value is also increasing when the effective attractive interaction 
J is decreasing. This situation is, in principle, similar to the case with attracting electrons, 
when there is no antiferromagnetic background for the carrier motion (see, for example 
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FIG. 1: The gap as function of carrier density at £3 = 1, and different values of t% and J. Here 
and below all functions are expressed in units of £3. Here and in Fig. 2b), we don't present some 
results in the case when the gap or the critical temperature is very small, since it was difficult to 
get accurate numerical solutions in these cases. 

IV. TEMPERATURE-CARRIER DENSITY PHASE DIAGRAM 



In this Section, we study the finite temperature properties of the model. It is known 
that there are two critical temperatures in the two-dimensional superconducting systems: 
the mean- field critical temperature T C MF , below which the uncorrelated pairs start to form, 
and the Berezinskii-Kosterlitz-Thouless critical temperature Tbkt < T^ F , below which 
the phases of the pair wave functions become algebraically ordered (for over-review, see for 
example 2lj]). Since such an order is the only possible order in our system, as it was stated 
above, we have to put T BKT = T c . 
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A. The mean-field critical temperature 



To find the dependence of the critical temperature T C MF on the particle density in the 
d-wave pairing channel, we need to solve the finite-temperature version of the set of Eqs. (|25j) 
and (|HJ). These equations follow from Q2D , © and (1221): 



'£ 2 (q) + AhJ(q) 



1 + tanh 



v / ^ 2 (q) + ^7, 2 (q)' 



e(k) 



(28) 



(29) 



2T J y/e*(k) + XiUM 

The system of equations for T C MF and // can be obtained from Eqs. (|2Hj) and (|29*j) by 
putting the amplitude of the order parameter equal to zero. The solution of these equation 
shows that the doping-dependence of the mean-field critical temperature (Figs. 2-4) have 



qualitatively ^ 



£ „ „ _ 




FIG. 2: The mean-field (left) and the BKT (right) critical temperatures as functions of carrier 
density at t2 = 1, £3 = and different values of J. 




FIG. 3: The mean-field critical temperature as a function of carrier density at £3 = 0.5. 
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FIG. 4: The same as in Fig.3 at t 3 = 1. 

It is important that superconductivity is suppressed at low carrier densities in the <i-wave 
pairing channel. To understand this qualitatively, one can notice that at low doping, the 
function 7^(q) in the equation for the mean- field critical temperature f)28|) can be approx- 
imated as 7 q ~ q 2 cos(2y? q )/2 ~ m*eF cos(2y9 qF )/2, where €f is the Fermi energy and m* 
is an effective quasiparticle mass, which is directly connected with i 2 an d t s (see the next 
Subsection). In this case, Eq. ()27|) has the form of a standard BSC equation for the critical 
temperature in the s-wave pairing channel with the renormalized coupling J — ► Jm* 2 e F /4. 
Therefore, the mean field critical temperature is exponentially small function of square of 
the carrier density S 2 ~ ep in this case. 



B. The critical temperature 

In order to study the BKT transition in the system, it is necessary to consider the 
superconducting order parameter transfer phase fluctuations. The phase of the order pa- 
rameter can be introduced in analogy with the fermion case, where the phase 6 nm of the 
fermion operator can be presented as a product of the neutral operator x{ T i n ) an d its phase 
exp{i#(r, n)/2} [21]. In our case: 

X 2 /' 2 = (X 1 ^ = X (r, n) exp{^(r, n)/2}. 

In this case, the superconducting order parameter can be presented as a product of its 
amplitude and the phase: 

(XWX* 1 ' 2 ) = 5 n , m+p A nm exp(z0 nm ). (30) 
It can be shown that the phase-dependence of the thermodynamic potential of the fermion 
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system with the Green function (}2*3"j) and the gap function (J3"U|) in the limit of small fluctu- 
ations of the phase of the order parameter is equal to 



n(A,0) = ^ J d 2 r{V6f 



(31) 



where the stiffness in the long-wave limit is 

5 1 1 f d 2 k 



4m* 16m* 2 T J (2tt) 2 cosh 2 [^( k ) 2 + A 2 7 2 (k)/2T] ' 



(32) 



22|). In our case, the effective mass of the free quasi-particles is m* 



3 -- 

(see, for example 
l/[4(t 2 + 2t 3 )]. 

In analogy with the 2D spin XF-model |23( , the equation for the BKT transition critical 
temperature, below which the phases of order parameter (the spin V6 orientation in the 
XF-model case) become algebraically ordered has the following form: 



7T 



T c = -J(A d ,fi',T c 



(33) 



where function J is defined in f)32p . 

The doping dependence of the superconducting critical temperature T c can be found by 
solving the system of Eqs. (|2*5|) . (|2*9Jl and (|3*3*j) . The solution of this set at different values 
of interaction shows that the doping-dependence of T c has qualitatively the same form as 
the doping-dependence of T 6 MF (Fig.4). It is possible to study some limiting cases of the 
solution of Eq. ()33|) analytically. In particular, in the limit of rather large carrier densities, 
when Arf <C T c (T c is close to T C MF ), one can make an expansion in powers of A^/T c on the 
right hand side of Eq. (J3*3~|) . In this case, this equation transforms to: 



T c 3 = AA 2 d (T c )5 3 /m* 



(34) 



where A ~ 7r 3 /128. Since the gap parameter depends on the critical temperature as A^(T) = 
Arf(0)[l — (T/T^ F ) 2 ] 1 ' a , where a > 1, at temperatures close to the mean-field critical 



temperature 



2J, the solution of Eq. is: 



T ~ 

- 1 c — 



T. 



MF 



1 

1 - - 

2 



*rpMFZ 



AA 2 (0) 



a/2 



pa/2 



(35) 



In other words, the critical temperature approaches to the mean-field critical temperature as 
the doping increases at large carrier densities. It is interesting, that equation is also valid 



13 



in the limit of low carrier concentrations when both critical temperatures are suppressed. 
Therefore, the second term in Eq. (J52j) is important at any carrier concentration in the d- 
wave pairing channel, contrary to the s-pairing case, where this term can be omitted at low 
values of 5 and T c ~ 5. 

It is important that the amplitude of T c is much smaller than T C MF , and the pseudogap 
region T c < T < T C MF is rather large in this case (Fig. 4). This is similar to the phase 
diagram of cuprates, where there is a large pseudogap region above the critical temperature 
at low carrier densities. It must be stressed, that in the case of higher carrier densities, one 
must condsider a model with more complicated effective interaction, comparing to Eq. (JBjl . 

In order to get superconductivity to develop at 5 ~ 0.05, similar to some of HTSCs, one 
needs to choose a pretty large interaction energy J > 5t%, which correspond to the energy 
scale of leV. However, to make a quantitative comparison with experiments, one needs to 
take into account different effects, which were omitted in this paper. In particular, optical 
phonons can give a significant contribution into the electron-electron interaction in cuprates, 
which can lead to decreasing of the required values of J in order to get superconductivity 
at 5 ~ 0.05. 



V. CONCLUSION 



To conclude, we have studied the superconducting properties of an effective model intro- 
duced in [n| in order to describe low carrier density properties of HTSCs. It was shown that 
the (i-wave pairing superconductivity in this model exists when the carrier density is larger 
then some critical density. This critical density strongly depends on the interaction energy, 
and it is growing with interaction decreasing. The amplitude of the pseudogap temperature 
is much larger than the critical temperature values, which resembles the experimental sit- 
uation on cuprates. There are some open issues which must be resolved. First of all, it is 
important to understand how to generalize the results on the case of larger carrier densities. 
In this case, the antiferromagnetic sublattice breaks down and it is not enough to put the 
attraction to be equal to the antiferromagnetic bond energy J. The effective antiferromag- 
netic attraction decreases. Also, the doping increasing is accompanied by increasing of the 
number of scattering centers created by dopants, which also leads to suppression of super- 
conductivity (see, for example Another important problem which is widely discussed 
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nowadays is to understand, whether the model can give an inhomogeneous superconducting 
state (see, e.g. |5] and references therein). It is also necessary to estimate pairings in other 
channels with different symmetry of the order parameter. These and some other questions 
are planned to be studied in the nearest future. 
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